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Abstract: We construct and study a dimensionally reduced effective tlieory for high- 
temperature SU(2) Yang-Mihs theory that respects all the symmetries of the underlying 
theory. Our main motivation is to study whether the correct treatment of the center 
symmetry can help extend the applicability of the dimensional reduction procedure towards 
the confinement transition. After performing perturbative matching to the full theory 
at asymptotically high temperatures, we map the phase diagram of the effective theory 
using non-perturbative lattice simulations. We find that at lower temperature the theory 
undergoes a second order confining phase transition, in complete analogy with the full 
theory, which is a direct consequence of having incorporated the center symmetry. 
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1. Introduction 

In the study of high temperature gauge theory in thermal equilibrium, a particularly useful 
approach has turned out to be that of dimensional reduction [1]. There, one describes the 
system via a d — 1 dimensional static effective theory built using the fact that, at high 
temperature, the non-static field modes decouple quite efficiently from the dynamics of 
length scales larger than or equal to the inverse Debye mass. This has led to several 
important advances, including e.g. new results for the spatial string tension and various 
correlation lengths [2], as well an efficient reorganization of the weak coupling expansion 
of the QCD pressure [3]. For the latter quantity, dimensional reduction has even provided 
a framework for extending the expansion to the full order, which is where the pressure 
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obtains its first non-perturbative contributions — a topic that has recently received much 
attention [4]. 

Despite the plethora of results derived utilizing dimensional reduction in perturbative 
setups, there is as such no reason to view the method as only a tool for weak coupling 
expansions. After all, it is well-known that at least a modest separation between the hard 
(27rT) and the soft and ultrasoft scales {m,i) , m^a_g) exists all the way down to the phase 
transition region. A minimal extension of the current perturbative approaches would be a 
non-perturbative, numerical determination of the couplings of the dimensionally reduced 
theory as a function of temperature, as attempted in Ref. [5]. Another more radical — 
and also more promising — approach is to attempt to cure the lack of center symmetry 
within the usual three-dimensional effective theory, EQCD, as it is the restoration of this 
symmetry that drives the phase transition to the confining phase in pure SU(A'') Yang-Mills 
theory. 

Recently, a new effective theory for high temperature SU(3) Yang-Mills theory has 
been proposed that has the virtue of respecting the Z(3) center symmetry of the original 
theory [6]. By construction, it cures the unphysical properties of EQCD that invalidate 
its use close to the dcconfinement transition region, such as the fact that the physically 
relevant phase of EQCD is in fact only metastable [7]. The new theory is thus expected 
to be valid over a wider temperature range than EQCD, and its properties, including the 
detailed structure of its phase diagram, are currently under non-perturbative study. The 
lattice simulations of the new theory are, however, technically quite demanding due to the 
high dimensionality of the respective parameter space [8] . 

To test the general idea behind the construction of dimensionally reduced center sym- 
metric theories, our aim in this paper is to perform an analysis of high temperature SU(2) 
Yang-Mills theory and its corresponding Z(2) invariant effective theory. The physical setup 
here is highly analogous to that of the SU(3) case, but has important simplifications due 
to the less complicated structure of the gauge group. Instead of having to introduce ten 
new heavy (and strictly speaking unphysical) degrees of freedom to ensure that a super- 
renormalizable theory with the correct symmetries can be constructed, it will this time be 
sufficient to consider only one extra real scalar field, whose mass is then the only adjustable 
parameter in the theory. This makes a marked difference from the point of view of the 
lattice simulations. 

A competing — or rather complementary — approach to the building of center- 
symmetric effective theories for thermal SU(A/") gauge theory is to set the requirement 
of renormalizability aside, but rather concentrate on the various kinds of models one can 
construct for the full theory Wilson line near Tc [9, 10]. An important step in this direction 
was taken in Ref. [11], where a new model was introduced that in a rough sense corre- 
sponds to a non-linear sigma model version of that of Ref. [6] . Very recently, a numerical 
study of this model was completed in the case of the gauge group SU(2) [12], and it will 
be interesting to investigate, to what extent their predictions will agree with those of our 
new theory. This task is, however, outside the scope of the present paper. 

The paper is organized as follows. In Section 2, we introduce the new effective theory, 
write down its Lagrangian and discuss its degrees of freedom as well as the structure of 
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its potential. We also analyze the vacua of the theory by computing the one-loop effective 
potential. In Section 3, we then perform a perturbative matching of the effective theory 
parameters to the full theory ones. This amounts to demanding that upon integrating 
out the extra heavy degree of freedom of our theory in the high-temperature limit, we 
obtain the usual EQCD Lagrangian, as well as making sure that the leading order tension 
and profile of the Z(2) domain wall in the effective theory agree with the corresponding 
quantities in the full theory. Section 4 on the other hand contains the details and results 
from a non-perturbative analysis of the new theory. We map its phase diagram in terms 
of the two remaining parameters, the four-dimensional gauge coupling and the ratio of the 
mass of the heavy scalar field and the temperature, and propose a simple scheme to match 
the value of the latter to the full theory. Finally, we summarize our findings and discuss 
their implications on future work within the present theory and its SU(3) counterpart in 
Section 5. Many computational details, such as the lattice-continuum relations for the 
parameters of the new theory, are presented in the Appendices. 

2. The effective theory 

Our aim in this Section is to generalize the center-symmetric effective field theory intro- 
duced for SU(3) Yang- Mills theory in Ref. [6] to the case of the gauge group SU(2). We 
begin this by introducing the required degrees of freedom and the Lagrangian, and then 
proceed to evaluate the effective potential of the theory to one loop order. Finally, by 
rewriting the Lagrangian of the theory in terms of fields fluctuating around the minima of 
the potential, wc make the connection to EQCD transparent. Our approach is somewhat 
different from that of Ref. [6], in which the effective potential turned out too complicated 
to evaluate. As will be seen below, having an explicit expression for Ves at hand makes 
in particular the perturbative matching of the theory to the full one intuitive and easily 
tractable. 

2.1 Degrees of freedom 

The degrees of freedom in any center symmetric effective theory, with which one wishes 
to correctly describe the physics of Yang-Mills theory close to the phase transition region, 
must include something reminiscent of the Wilson line of the full theory. 



Denoting this 2x2 matrix field in our theory hy 2, we require that the action of the theory 
must be invariant under the local gauge transformations 



with A(x) denoting the spatial gluons and s(x) G SU(2), as well as under the global Z(2) 
phase rotations. 




(2.1) 



2:(x) ^ s(x) Z{:x.) s(x)t , 

A(x) ^ s(x) (A(x) + iV) s(x)t , 



(2.2) 




e'''".E(x) = ±Z(x) 



(2.3) 
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for integer n. 

In analogy with Ref. [6], we do not wish to formulate our theory in terms of an SU(2) 
matrix, but instead regard the field Z as a spatially coarse-grained Wilson line operator, in 
terms of which a super-renormalizable effective theory with only polynomial interactions 
can be constructed. This implies defining Z via the block transformation 

Z{^) = -^ f d3yi7(x,y)0(y)C/(y,x), (2.4) 
v'Block Jv 

where the integration goes over the (somewhat arbitrary) 0{T''^) volume of the block and 
L'"(x, y) is a Wilson line connecting the points x and y at constant time r = 0. 

The crucial difference to the SU(3) case now comes from the fact that with SU(2), the 
coarse-graining almost preserves the group property, as a sum of SU(2) matrices is up to 
an overall real constant another SU(2) matrix. Taking furthermore into account that the 
exponentiation of a sum of Pauli sigma matrices, or the generators of SU(2), can be written 
as a linear combination of the very same matrices plus the unit matrix, this implies that 
we may parametrize Z as 

z = ^^^m + iUaaa} ^ + m. (2.5) 

Here, S and Ha, a= 1, 2, 3, are real scalar fields, the aa denote the Pauli matrices, and the 
new fields are obviously given by the formulas 

E = TrZ, (2.6) 
m = Z-^TTZl. (2.7) 

Compared to the Wilson line operator of the full theory, the field Z seems to contain one 
extra real degree of freedom. In analogy with Ref. [6], we however assume this auxiliary field 
to be heavy (with a mass of order T that originates from the coarse-graining procedure), 
so that at high temperatures it decouples from the physics of length scales > 1/ {gT) and 
can be integrated out. We will identify this heavy mode later. 

2.2 The Lagrangian 

We take the Lagrangian of our effective theory to be composed of the conventional kinetic 
terms for the magnetic gluons and the Z field, plus a potential V{Z), 

£ = 53-2 { i TV 4 + IV ( A-Z^ A-Z) +V{Z)Y (2.8) 

with Di = di — i[Ai, ■ ] and Fij = diAj — djAi — [Ai, Aj]. The potential V is chosen to be 
of the most general superrenormalizable, Z{2) and gauge invariant form^ 

V{Z) = + b2Ul + ciE^ + C2 {Ulf + caE^n^, (2.9) 



^We want to emphasize here that unhke in the SU(3) case, our potential contains all terms allowed by 
symmetry and supcr-rcnormalizability, so that this time even the purists should not object to calling the 
theory a true effective theory rather than just a model. 
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where, due to the superrenormahzabihty of the theory, the coefficients Cj do not depend on 
the renormahzation scale of the effective theory. The mass terms 61 and 62, on the other 
hand, acquire a scale dependence at two-loop level. In the MS scheme, their expressions 
read 

bi{fl3d) = [48c? + 124 + Ucsgl] log (^^^ , (2.10) 

b2{fl3d) = [80ci + 4ci - 40c25i] log (|^) > (2-11) 

where /xsd is the renormahzation scale of the effective theory and A^. are constants speci- 
fying the theory. In the following, we will set fisd = 93 and abbreviate 61(53) = bi- 
We can write the potential also in the alternative form 

V{Z) = hiTr{Z^Z) + h2iTrZ^Zf + gl\^siTiU^ + S2{TiU^)^ + S3^^], (2.12) 

where we, motivated by a perturbative analysis, have assumed the coupling constants hi 

and Si to scale as (gs)^, so that the operators become naturally divided into hard and soft 
subsets. The reason for this is related to the transformation properties of the operators 
under an SU(2)xSU(2) transformation of Z, 

z nizn2, ^i e su(2), (2.13) 

which will later be seen to translate into a shift symmetry in the light physical field of 
the theory, when we integrate out the auxiliary heavy scalar. The SU(2)xSU(2) invariant 
terms in the potential are responsible for the 0{T) mass scale of this scalar field, but do 
not affect the perturbative matching to EQCD at leading order. They must therefore come 
with sufficiently hard, or order {93)^, couplings denoted by hi. The part of the potential 
that does not possess this additional symmetry will, on the other hand, be found to give 
rise to the potential and mass scale of the Goldstone mode of the hard potential, later to 
be identified as the Aq field of EQCD. Consequently, these terms in the potential come 
with 0{g^) couplings, parametrized by the constants Sj. Finally, relating the two sets of 
couplings to each other produces 

61 = ^^1, b2 = ^{hi+glsi), (2.14) 

Cl = ^^2+53^3, C2 = ^(/l2+5is2), C3 = ^/l2. (2.15) 

Perturbatively matched results for these parameters will be listed in Section 3.3. 



2.3 The effective potential 

Before moving on to the perturbative matching of the effective theory to the full one, 
we address a somewhat more general question that will turn out to be of much use in the 
forthcoming analysis: What is the one-loop effective potential of our theory when computed 
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in the background of constant or sufficiently slowly varying^ classical fields 

(S) = p, (2.16) 

IJia) = Oj6a,3 . (2.17) 

The choice to have a non-zero expectation value only for the a = 3 component of the 11 
field is totally arbitrary, but naturally allowed by gauge invariance. Likewise, we could 
have included a non-zero {A'^), but left it out for brevity, as it will not be needed in the 
present context. 

A straightforward calculation (the details of which are given in Appendix A) gives for 
the effective potential in the gauge 



VeS = gf {hp^ + b20J^ + Cip^ + C2UJ^ + C^pW) - ^2 - f^^ 

^{{bi + b2 + Gp^ci + 6u^C2 + {p^ + w2)c3 - ^/v) 



Gvr 

1 (;i -I- no -t- K)u (;i -h \)LU co -t- \ u ^ uj iii-i — a / // 

127r 

+ (61 +b2+ ep^ci + 6io'^C2 + {p^ + w2)c3 + 



+ 2 (ec^2 + 262 + 4u;2c2 + 2p'c^f^ } + 0{gl), (2.18) 

where we have denoted 

?7 = (61 - 62 + Q{p^ci - uj^C2) - {p^ - uj'^y^f + IQp^oj^cl (2.19) 

The explicit gauge parameter dependence of the result should come as no surprise, as the 
effective potential is not a gauge invariant quantity. Only its minima structure is of physical 
significance. 

2.4 Minimizing the potential 

The truncation of the loop expansion of the effective potential to one loop order implicitly 
assumes that we are working in the weak coupling limit, so we may consistently further 
simplify the result of Eq. ( p. 18 ) by using the identities of Eqs. (2.14)-(2.15). Looking at 
the potential order by order in (73, we first obtain 

Kfr = ^ (P' + '^') (2/ii + h2{p' + ^')) +0[{g^f) , (2.20) 

which, depending on the sign of hi (note that perturbative stability requires /i2 > 0), is 
either minimized by p = a; = or 

p' + u^' = v' = -^+0{gl). (2.21) 
112 

Upon comparison with the effective potential of the Wilson line in the four-dimensional 
theory, we note that the case of interest in the present paper, where we wish to be working 



^Implying that we may neglect the derivative terms of the background fields in the classical action when 
computing the one-loop functional determinants in Appendix A. 
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in the high temperature (deconfined) phase of the full theory, is that of hi < 0. We will 
thus assume this to be the case from now on. Using then the fact that any 0{g^) correction 
to the functions p and uj would only contribute to the potential at the NNLO level, we may 
parametrize p and oj by 

p = vcos{7ra), (2.22) 
Lo = vsm{Tra), (2.23) 

with a a real function, and insert these values into Eq. (|2.18|). This leads to 



2 4 3 

Ves = ^ sin2(7ra) + ^ sin^(7ra) + ssv^ cos^{-Ka) - ^ \ sin(7ra)|3 + 0{g1), (2.24) 
2 4 Svr 

where we have dropped a trivial constant from the result. Note that the ^-dependence of 
the effective potential T4ff in Eq. ( |2.18| ) has disappeared at the minimum of the potential, 
defined by Eq. ( 2.21 ), as required by Nielsen's theorem. 

The locations of the minima of the function V^s in Eq. ( |2.24 ) depend on the values of v 
and Si. As we want our theory to inherit the Z(2) minima structure of the effective potential 
of the full theory Wilson line, we will for the time being assume (and later confirm) that 
these parameters have values such that the potential is minimized by a = rnr, n G Z. This 
in turn implies that the effective potential has its minima at a; = 0, or 

{Z) = ±^l. (2.25) 

Specializing now to fluctuations around one of these physically equivalent Z(2) minima, 
we denote 

^vl + gsl^^^l + ix)^ (2.26) 

with X 3- traceless hermitian matrix field, x = Xa^^af^, and in addition rescale the three- 
dimensional gauge field by Ai — > g^Ai. This enables us to rewrite the Lagrangian of the 
effective theory in the form 

£= ^TrF^j + ^[{di<t>f + ml^^]+TT[{D,x? + mlx']+Vs{c^,x), (2-27) 
where the masses of the shifted fields have become 

ml = ml = Sv'^ci = -2hi, (2.28) 
ml = 2 (62 + v^c^) = {si - 4i;2s3) gl (2.29) 

The parametrically heavier field (p is hereby identified as the auxiliary degree of freedom 
in our theory, which is to be integrated out in the high temperature limit. Finally, given 
in terms of the shifted fields, the potential Vg reads 

Vs = 2g^v[2ci(t)^ + C3<AXa} + 5i{ci<A' + C2 [xlf + c^4>\l\ (2-30) 

= g,vh,4>{ct>' + xi) + '^ (<A^ + xif + ^9lvs,<f' + f {^s,4>' + .2 {xlf ]■ 
The Z(2) symmetry is now realized by a simultaneous sign flip of v and (p. 
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3. Parameter matching at high temperature 



In the previous Section, we saw that at least under certain assumptions on the values of 
its parameters, our effective theory possesses a minima structure reminiscent of that of 
the underlying full Yang-Mills theory. Motivated by the connection between our degree of 
freedom Z and the Wilson line operator Q, it is therefore tempting to associate our light 
field X) describing traceless anti-Hermitian fluctuations around the deconfining minima, 
with the Aq field of EQCD. To make this correspondence more concrete, we now move on 
to the high-temperature limit, where we integrate the heavy fluctuation field out from the 
effective theory and then match the Lagrangian of the resulting 3d Yang-Mills plus adjoint 
Higgs theory to that of EQCD. Following the approach of Ref. [6], we will also require that 
the leading order Z(2) domain wall stretching between the two deconfined minima of the 
effective theory has the correct tension and width dictated by the full theory. It will be 
seen that this is sufficient to fix all but one of the coupling constants of our theory. 



3.1 Reduction to EQCD 

The process of integrating out the heavy field (/> is greatly simplified by the SU(2)xSU(2) 
invariance of the hard part of our potential, as it is easily seen to translate into a shift 
symmetry in x (for details, see Ref. [6]). This implies that upon integrating the heavy field 
out, the graphs with only /i-vertices will have to cancel, and, in particular, that at leading 
order we may simply read off the parameters of the resulting Lagrangian, 

Aight = \ Tr 4 + Tr[(A xf + + + -, (3.1) 

from Eq. ( 2.27] ). Doing so, we immediately obtain 



ml = {si - 4vhs) gl + O [gt) , ~X = 2s2 gl + o(^) , (3.2) 

where we have used the SU(2) identity (Trx'^)^ = 2Trx'*. A comparison with the EQCD 
Lagrangian then confirms the validity of the simple identification Xa ^ further 
leads to the leading order parameter values 

,2 



mo 



93 + 0[^) =g'T + 0{g^T), (3.3) 
si - 4vhs + O {gl) = ^ + 0[g^T) , (3.4) 



2T 
"3 



where g is the four-dimensional gauge coupling and T the temperature of the full theory. 
It is straightforward to see that these results are consistent with the assumption we made 
in Sec. 2.4 about the 0(y{g^)^) effective potential being minimized by vra. 
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Figure 1: Left: The normalized domain wall free energy densities T{z) in the effective and full 
theories, plotted versus z = g{T)Tz. Right: The relative difference of the two curves, {Tcs{z) — 

3.2 Domain wall properties 

We now move on to the determination of the leading order Z(2) domain wall profile in the 
effective theory, stretching between the two physically equivalent Z(2) minima of Eq. ( 2.25| ). 



It is clearly possible to have the field Z minimize the hard part of the potential throughout 
the wall, which implies that the width of the wall will become of order \/{gT) and, in 
particular, that we may use the results of Sec. 2.4 for the perturbative effective potential of 
our theory.^ In Appendix B, we present in some detail our leading-order calculation of the 
domain wall profile, including the gauge- invariant free energy density J'{z), from which we 
can extract the domain wall tension (by integrating J-{z) over z) and width (by integrating 
z'^!F{z)). These two quantities agree with those determined perturbatively from the full 
theory at high temperature, if one chooses the parameters as 

V K, 2.001622T, (3.6) 
S3 K, 1. 60379/ (127r2r). (3.7) 

We observe from Fig. 1 that with this choice, the wall profiles in the two theories are 
practically overlapping, with the relative error being in the per cent range. 

3.3 Leading order parameters 

To summarize our findings from the present Section, we collect here all the results from 
the leading order perturbative matching of the effective theory parameters. Defining the 
dimensionless ratio thq/T = = r, with v = v/T, we may write the couplings in 



The assumption we made in the derivation of Eq. (2.18) about the background field being sufficiently 
slowly varying in order for its kinetic terms to be negligible is found to be fulfilled here, as the width of 
the wall being of order l/(gT) implies that each spatial derivative of the background field introduces an 
additional factor of g. 



terms of only g, r and T, giving 

1 



61 = -^r^T^ (3.8) 

62 = -^r^T^ + 0.44184l52r2, (3.9) 

ci = 0.0311994r2 + 0.01354155^ (3.10) 

C2 = 0.0311994r2 + 0.008443432^2, (3.11) 

C3 = 0.0623987r2, (3.12) 



and somewhat trivially 



gl = g'T. (3.13) 
These results will be put to use in the following Sections. 

4. Beyond perturbation theory 

Having the perturbatively matched theory at hand, an immediate question to ask becomes, 
what happens to its properties when one leaves the weakly coupled regime we have been 
investigating so far. To answer this, we now move on to consider non-perturbative numerical 
simulations of the theory defined by Eqs. (^), (2.9) and (^)- (|3.13| ). Ideally, one would 



of course like to be able to analyze the theory in all generality, i.e. in terms of completely 
general couplings 6j, q, but as we will see in the following, the two-dimensional slice spanned 
by r and g already provides a physically motivated, interesting and, above all, numerically 
controllable subset of the entire parameter space. 

4.1 Numerical simulations 

In the present paper, we wish to address the physically perhaps most pertinent question of 
what the phase structure of our theory looks like in the (r, g) plane. Prom the evaluation 
of the one-loop effective potential, we already know that for any finite real value of r, the 
theory will always find itself in the deconfined phase in the limit of 5 — > 0, but in order to 
look at non-zero values of 5, we must formulate it on a discrete space-time lattice. This 
procedure consists essentially of two parts: First, writing down the action of the theory 
in a discrete form such that it reproduces the same long distance physics that one would 
obtain in the continuum MS renormalization scheme and, second, performing numerical 
simulations. 

The first goal can be achieved by naively discretizing the action and then matching 
the parameters of the discretized theory to the MS ones, imposing the condition that the 
pole masses and interaction vertices of the two theories coincide. The matching can be 
performed to any desired order in the lattice spacing a, due to the superrenormalizability 
of the theory, within the framework of lattice perturbation theory. Here, we perform the 
calculation up to 0{a) (or equivalently 0(l//3), where /3 = A/{ag'^) is the lattice coupling), 
thus ensuring that all non-perturbative properties of the lattice theory coincide with those 
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Figure 2: The expectation value and susceptibility — as fmictions of at 

= 5 with (/3, A^) = (12, 64) [left] and (6, 64) [right]. The critical region shrinks as the physical 
volume is enlarged. 

of our continuum effective theory in the continuum limit /3 ^ oo. For further technical 
details, we refer the reader to Appendix C and the corresponding calculation performed 
in the SU(3) case in Ref. [8]. Our lattice action is defined by Eqs. ( C.1 )-( C.8 ). Note the 
manifest Z(2) symmetry S <-> — S. 

In our simulations, we have used the Metropolis algorithm for updating the S and 
n fields, and the Kennedy-Pendleton heatbath and overrelaxation algorithms for the link 
variables. At small enough g, the theory is seen to reside in the deconfined phase, but for 
each r, there is some non-zero value of g, at which the system exhibits a Z(2)-restoring 
phase transition to the confined phase, characterized by the vanishing of the expectation 
value (S). There is no visible discontinuity in (S), which suggests a second order transition, 
as expected from universality arguments (Z(2) symmetry at d = 3). For details, see Fig. |2|. 

Next, we determine the pseudo-critical point gdr) from the peak of the susceptibility, 

X. ^ m')-m)\ (4.1) 

fitting the three highest points with a second-order polynomial and estimating the error 
using ten jackknife blocks. In the thermodynamical limit, gdr) is known to coincide with 
the physical critical point. To further ascertain the universality class of the transition, we 
measure the Binder cumulant = {T,^)/{T,^f. In Fig. |, we exhibit our measurements 
for several lattice sizes as functions of the rescaled variable (l/g"^ — l/g'^)N^^'^ , where we 
take the d = 3 Ising value 0.63 for the critical exponent ly [13]. A satisfactory collapse of 
the data is observed at = 5 and 10 (see Fig. |3|), and the value of the cumulant aX g = gc 
is consistent with the Ising value 1.604. 

To obtain the phase diagram of the continuum theory, we must perform thermodynam- 
ical and continuum extrapolations, Na — > oo and a — > 0. The leading finite-size corrections 
to the critical coupling g^. are described by standard renormalization group analysis, and 
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Figure 3: A check of the universaUty class of the transition: The Binder cumulant 
is shown as a function of the rescaled variable {1/ — Xj g^^N^I" for — 5 (left) and 10 (right). 
The pseudocritical coupling is the value of the coupling, which maximizes the quantity Xs in 
Eq. (41), and v = 0.63 as appropriate for a three dimensional Z(2)-transition. A satisfactory data 
collapse is observed for various volumes, and the cumulant value at is consistent with the 3d 
Z(2) value 1.604. 



are of the form 

l/g^Na) = 1/52(00) + C(iVa)-V-, (4.2) 
where gc{Na) is the pseudo-critical value of the coupling which maximizes the susceptibility, 



Eq. (4.1), on a lattice of size Na, and C is a constant fitted to our data. This ansatz 
describes our large-volume data well, as illustrated in Fig. ^. Note, however, that the 
finite-size effects increase as r is reduced, which is no surprise, as r determines the mass 
mo of our heavy degree of freedom, whose Compton wavelength increases as 1/r. This 
means that very large lattices are required to study the regime r ^ 1. 

Turning then to the continuum extrapolation, Fig. |5| shows our measurements of 1/(7^ 
in a small, constant physical volume {Na{l3)Y as functions of 1//3. Our data are consistent 
with a discretization error linear in 1//3, i.e. proportional to the lattice spacing a. The 
origin of this 0{a) error can be traced back to the imperfect matching of the parameters 
of the lattice potential with their continuum counterparts, displayed in Eqs. ( p.l| ): An 
0{a) error remains, which would take a three-loop calculation to eliminate. Therefore, we 
extrapolate our data linearly to 1//3 = 0. In any case, the discretization error is small, 
which we expect to remain true as long as all physical length scales of our theory are 
sufficiently large compared to the lattice spacing: 

a^{l/{g'T)M{gT)M{rT)}. (4.3) 

The last of these inequalities limits our ability to study the regime of large r. Therefore, we 
have concentrated our simulation efforts on = 5 and 10, for which the thermodynamical 
and continuum extrapolations are both under good control. 
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Figure 4: The behavior of the pseudo-critical couphng \/g1 for fixed = 5 and 10, as one 
increases the extent of the lattice, N . The solid and dotted curves are the renormalization group 
scaling fits to the four and three rightmost datapoints, respectively, and the two dashed horizontal 
lines correspond to the iV — > oo limits of these curves. The effect of including or neglecting the 
fourth data point is of the same order of magnitude as the statistical errors in the fits. We use this 
to quantify the systematics of our fits, and add the two sources of error in the final error estimate 
of our extrapolation. 

4.2 Phase diagram and the matching of r 

Having taken into account the errors caused by the finite values of the lattice size and 
spacing as explained above, we obtain the phase diagram of the continuum effective theory 
exhibited in Fig. ^. In addition to our two data points, we show the special point at r = 0, 
as obtained from Ref. [14]. While for r = our numerical accuracy is severely hampered 
by problems in the thermodynamical extrapolation, such simulations are not needed, as 
for r = the field S completely decouples from the rest. The relevant part of the action 
then becomes that of a three-dimensional one-component theory, the phase diagram of 
which has already been carefully determined in Ref. [14]. This reference provides us with 
an accurate value for the end point of the critical line, as explained below. 
Quoting Ref. [14], the critical mass of a theory defined by the Lagrangian 




(4.4) 



reads 



= A/3) 
A2 



0.0015249(48). 



(4.5) 



A conversion of this result to our theory and the renormalization scale fJ- = = g 



'^T then 
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Figure 5: The effects of the finite lattice spacing on the value of the pseudo-critical coupling 1/ gl 
in a smah fixed physical volume V = (16/3)^ /gf , with N) = (6, 8), (9, 12), (12, 16), and (18, 24). 
The efl^ect of finite lattice spacing is estimated from the slope of a straight fit to the data points 
giving l/gl = 0.287(7)+0.14(5)//3, with xVdof = 1.03/2, for = 5 and l/gl = 0.37(l)+0.21(8)//3, 
with xVdof = 0.96/2, for = 10. 



reveals that at r = 0, the critical value of our gauge coupling l/g^ reads'* 

l/c/clr=o = 0.025543. (4.6) 

This point fixes the form of our critical line at small r, and shows that there is a critical 
coupling of 5 ~ 6.3, beyond which the effective theory is always in the confined phase, no 
matter how small one dials the mass of the heavy scalar field. 

Moving on to non-zero values of r, we observe that the confined and deconfined phases 
in our theory are separated by a second order phase transition line, which goes through 
the data points (r^, l/g"^) = (5, 0.175(4)) and (10, 0.206(7)). Fitting a polynomial trial 
function to these points as well as to the known intercept of the transition line at r = 0, we 
obtain the curve shown in Fig. ^. One observes that the curve is rather flat for r ~ 0(1) 
and higher. This can be understood by noting that the phase transition is a long distance 
property of the theory, while the scale r originates from the microscopic coarse-graining 
of the Wilson line. Thus, while r is formally an undetermined parameter of our theory, 
we in fact expect (and here observe) that the properties we are interested in are rather 
insensitive to it. 

In particular, the latter statement implies that as r — > oo, the critical curve will 
approach a horizontal line, as the radial fluctuations of the field Z become frozen and 
Z simply reads v/2 times an SU(2) matrix. In this limit, the mass of the heavy field 
of our effective theory becomes so high that it decouples from the physics and can no 
longer affect the properties of the system, such as the critical value of However, 
the potential V{Z) in Eq. ( 2.12| ) is not 0(4)-symmetric. It has two degenerate minima at 



S = zbu, n = 0, separated by barriers of height 0(5(3). Thus, a Z(2) transition can still take 
place as a function of g^, leading to a deconfined phase where S = ±v at high temperature. 



*A consistent value = A/m ~ 24 has also been obtained in Ref. [15]. 
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Figure 6: The phase diagram of the perturbatively matched theory on the (r^, 1/g^) plane. The 
solid blue data points are from our numerical simulations and the open red one has been obtained 
from the known location of the critical point of A(/)^ theory, Eq. ( |4.6| ). The curve connecting the 
points has been obtained from a polynomial fit and has been included in the diagram to guide the 
eye. 



Finally, we would like to choose a value for the heavy mass parameter r, so as to best 
reproduce the properties of the original (3 + V)d SU(2) Yang-Mills theory. Ideally, one 
would like to perform a non-perturbative determination of all the couplings of the effective 
action, including r, as a function of the temperature of the (3 + l)d theory, by measuring 
the large-distance two-point functions of the relevant fields. As a first simple attempt in 
this direction, we will now retain our perturbative determination of the other coefficients 
and fix r by demanding that the transition temperature be the same in the full and the 
effective theories. Assuming further a one-loop running of the full theory gauge coupling 

with /2 = 6.74r [16], and using the result Tc/Aqcd ~ 1-23 [17], this amounts to fixing the 
coupling at the transition temperature to be g^i^c) ~ 5.1. From the phase transition 
line of Fig. ^, we then obtain 

r = 2.6(5), (4.8) 
which in this rather crude approximation scheme is our final result for the parameter. 
4.3 Discussion 

We have formulated our effective theory on a discrete three-dimensional lattice, which 
has enabled us to investigate its properties at finite couplings. The parameters of the 
theory used in the simulations were partially fixed by perturbative matching to the full 
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theory, which left only the parameter r, related to the heavy mass scale of the effective 
theory, unfixed. It is an impressive demonstration of the power of respecting the correct 
symmetries that starting from such restrictive assumptions, clearly not fully consistent 
close to the phase transition region, we were able to find a line of second order phase 
transitions in the (r, g) plane, just as the full theory predicts. 

The success of the new effective theory provides a stark contrast to the case of EQCD, 
which cannot reproduce the phase structure of the full theory, as it inherently only describes 
small fluctuations of the Aq field around one of the Z(2) minima. The two deconfining 
minima of this theory are therefore physically inequivalent and only resemble each other in 
the immediate vicinity of the critical point, which furthermore is located far away from the 
physical line of matching [18]. In addition, in the absence of an order parameter, there is no 
clear signature of a confining phase, whereas our theory provides a consistent description 
of it. 

Having said this, it should nevertheless be emphasized that there is still work to be 
done within the new effective theory. First of all, while the perturbative matching we 
have performed in the present paper ensures that the theory has the correct behavior at 
asymptotically high temperatures, one would like to quantify the extent to which it is able 
to capture the dynamics of the phase transition. To this end, one should compute various 
physical quantities in the effective theory and compare them to results in the full theory 
as a function of the temperature; among these are e.g. the spatial string tension and the 
domain wall tension, as well as the entire spectrum of the effective theory, which is related 
to the screening masses of the full theory. After the matching of r in Eq. ( [4. is taken 
into account, our theory is in principle fully predictive, so all of these comparisons will be 
non-trivial tests of the validity of the assumptions and claims we have made. 

While our effective theory should already at the present stage give a good approximate 
description of high-temperature Yang-Mills theory down to significantly lower temperatures 
than EQCD, there are also various ways to improve the approximation further, such as 
by considering a non-perturbative determination of the effective theory couplings. The 
simplest strategy would be to keep the present perturbative matching of all coefficients, and 
let g and r evolve non-perturbatively with T. The gauge coupling g can be determined from 
simple measurements of the spatial string tension or the domain wall tension as already 
done in Ref. [20]. The parameter r, or equivalently the mass of the heavy field which 
represents the magnitude of a coarse-grained Wilson line, can on the other hand be also 
determined in the full theory, by performing a blocking transformation on the Wilson line 
and monitoring the correlator of its magnitude. These two non-perturbative matchings are 
straightforward to implement, and should improve the quantitative description the effective 
theory provides at temperatures close to Tc. 

5. Conclusions 

Building a dimensionally reduced effective theory for high-temperature SU(A^) Yang-Mills 
theory that respects all symmetries of the latter is an important task not only from the 
point of view of completeness and rigor, but also because there is genuine interest to 
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study the extent to which the failure of EQCD to respect the center symmetry explains 
its unphysical properties close to the phase transition region. To this end, the task of 
supplementing EQCD with just enough heavy degrees of freedom to enable the formulation 
of a superrenormalizable center-symmetric theory was tackled in Ref. [6], where a new 
effective theory for SU(3) Yang-Mills theory was built. An unfortunate practical problem, 
however, turned out to be the high dimensionality of the respective parameter space, as 
well as the fact that not all operators allowed by the symmetries could be included in the 
potential for practical purposes. Consequently, numerical studies of the theory have turned 
out to be very impractical [8]. 

As conventional wisdom states that 2 ^ 3, it is not surprising that SU(2) Yang-Mills 
theory provides a technically much simpler platform to investigate the same issues one 
encounters with SU(3) Yang-Mills theory. In the present paper, we have transcribed the 
ideas of Ref. [6] to the case of the simpler gauge group, and doing this have managed to 
include in the Lagrangian all superrenormalizable operators allowed by the symmetries. 
In addition, we have seen that the number of heavy scalar degrees of freedom one needs 
to introduce in the theory as well as the number of free parameters not fixed by simple 
pcrturbativc matching arc both one. This has enabled us to extend the matching to 
include the determination of the mass scale of our heavy field, using as matching criterion 
the simple requirement that the phase transition of the effective theory take place at the 
same temperature as that of the full one. More generally, we have been able to map the 
phase diagram of the new theory as a function of the mass parameter r and the four- 
dimensional gauge coupling g, revealing a second order phase transition line as predicted 
by the full theory. This required performing non-perturbative lattice simulations, as well 
as a two-loop lattice perturbation theory calculation. 

Our primary motivation in the present paper has been to address the question, whether 
the region of applicability of dimensionally reduced effective theories of hot Yang-Mills 
theory can be extended down to the phase transition region, provided one correctly accounts 
for all the symmetries of the full theory. Our findings regarding the phase structure of 
our new effective theory are highly encouraging in this respect: The spurious (^o) 7^ 
phase of EQCD disappears, and we observe instead a genuine Z(2) finite temperature 
transition, which can be adjusted to occur at the right temperature. Nevertheless, more 
work is required to put our approach on solid quantitative ground. In this process, the 
framework laid out in the present paper will be an invaluable asset, as the theory is now fully 
predictive and simply awaits further simulations to be performed. Beyond this, obvious 
further challenges include the generalization of our work to the case of SU(3), as well as 
the introduction of soft center-symmetry breaking operators in the effective theory with 
the purpose of modeling the presence of dynamical quarks in the full one. 
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A. Calculation of the effective potential 

The evaluation of the effective potential proceeds in the standard way of adding to the tree 
level potential of Eq. (2.9) the contributions from the one-loop fluctuations of the S, 11 
and gluon fields. For this purpose, we write our scalar fields in the form 

S = p + 53</', (A.l) 

Ila = t^(Ja,3 +93Xa, (A.2) 

where (p and x will be referred to as the fluctuation fields from now on. The one-loop 
functional determinants are computed in the most straightforward way by choosing to 
work in the renormalization, or gauges, defined by the gauge fixing function 

Ga = + ^93^ {SalX2 " <5„2Xl))- (A.3) 

This enables us to decouple the gluons from the x field at the level of the quadratic action, 
leaving as the only non-diagonal operator the mixed mass term of and X3- 
The general structure of our one-loop potential now becomes 

Ves = \Vd + VA + Fgh + V^^ + V^, (A.4) 

^3 



where V^i denotes the potential of Eq. (p.9|), evaluated with the classical fields, and the 
remaining terms are the contributions from the one-loop functional determinants of the 
quadratic gluon, ghost, (p and x actions, respectively. The second-to-last term corre- 
sponds to the contribution of the entangled cj) and X3 fields, while V^^^ denotes the one-loop 
determinant coming from xi and X2- 

A straightforward calculation produces for each of the above functions 



Vci = hip^ + 62^^ + ci/ + C2UJ^ + cap^w^ (A.5) 
1 
2 



Vk = \ Tr ln<| ( - 5i,d'' + (1 - l/i) didAs'^'' + u:'' {Sa,iSb,i + 5a,2<56,2 



+ (A.6) 



= -Trln{ -d^ + ioj^ {5a,i6k,i + 5a,2<5fe,2) } 

■e'\ (A.7) 



37r 
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V^^ = ^ TVln{ - 92 + 61 + 62 + 6(/>2ci + oj^C2) + {p^ + uj^)c3 - Vv} 
+ ^ Tr Inj - + 61 + 62 + 6(/>2ci + uj'^C2) + {p^ + io'^)c3 + ^} 
= - ( (foi + 62 + 6(/52ci + u;^C2)c2 + (p^^ + u;^)c3 - ^f'^ 



+ (61 + 62 + 6(p'ci + L^2c2)c2 + + c^2)^^ ^ ^) 3/2 ^ 

= Trln{ - + + 262 + 4w^C2 + 2^203} 

= -7^ (e^' + 262 + 4a;'c2 + 2p2c3)3/2 ^ ^^_g^ 

DVT ' 

where we have denoted 

r/ = (61 - 62 + ^{p^ci - cj^ca) - (p2 - a;2)c3)2 + IGp^w^^i. (A.IO) 



Adding up the various contributions now finally leads to Eq. ( ^.ISl) . 
B. Effective theory domain wall 

In this Appendix, we compute the leading order profile of a domain wall that stretches 
between the two physical minima of our theory, Eq. (|2.25| ), our purpose being to find out, 
for which values of the parameters v and S3T the effective theory reproduces the domain 
wall tension and width of the full theory. 

Taking the wall to lie in the (x, y) plane and using the same gauge and parametrization 
of Z as in Sec. 2.4, we have 

Z{z) = I diag(e*™("), e-*™(")) , (B.l) 

where a is this time required to satisfy the domain wall boundary conditions 

a(z = -oo) = 0, a{z = 00) = I. (B.2) 

The wall Lagrangian is composed of the usual kinetic term for Z plus the one-loop effective 
potential, giving 

/:waii = 93^'Tr [OiZ^DiZ) + Feff(^), (B.3) 

where V^fj can be read off from Eq. ( |2.24| ). Using the results from Section 3.1, Eqs. ( |3.3| )- 
( |3.5| ), and defining 

V,^{Z) ^ -^U{a), (B.4) 

we therefore see that the domain wall tension, i.e. its free energy per unit area, can be 
written in the form 

1 

^dw [a] = - 53 ' (^^ Tf t3/2 / dz{ia'f + U{a)}, (B.5) 

^ J-OD 
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with 

U{a) = sin^(7rQ;) + (l + 127r^rs3) sin^(7ra) - |sin(7ra)| (B.6) 



z = gTz , ^ = • (^•'^) 



and 



Minimizing this functional, subject to the boundary conditions of Eq. ( |B.2| ), will yield the 
leading order domain wall profile and tension. 

Following the steps outlined in Ref. [6], it is easy to find that the classical equation of 
motion for a{z) can be written in the form 



'1/2 

while the resulting domain wall tension is 



a = FdwH = {-KvTfT^/^ [ da U{a] 

Jo 



1/2 



(B.9) 



We define the width^ Az of the domain wall as the square root of the ratio of the second 
moment of the domain wall free energy density and the tension, which gives 

(Az)2 ^ r dzz' U{a{z)) 

93 J-oo 

^2{7wfT^ /" daiUia^fl^i H da^U^a^T''^^ . (B.IO) 

ffS^" Jl/2 \Jl/2 J 

The integrations in these three equations are trivial to perform numerically for different 
values of the unknown parameters v and S3T. 

In high temperature SU (2) Yang- Mills theory, the domain wall free energy density has 
been found to equal [19] 

resulting in the domain wall tension 
and width squared 

/CO ^2 /2 3 

dz z^ J'Yuiz) = ^^^^^ . (B.13) 



^Analogously to the Z(3) case [6], we have observed that different definitions of the domain wall width 
lead to almost identical values for v and S3. This can be attributed to the highly similar potentials (and 
subsequently similar wall profiles) in the effective and full theories, and may be taken as an a posteriori 
justification of leaving out the non-superrenormalizable terms from the effective theory Lagrangian. 
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A straightforward numerical calculation shows that the effective theory reproduces the 
latter two results with the parameter values given in Eqs. (|3.6D "(^). Using these values, 
the energy density profiles of the wall in the full and effective theories are found to agree 
within a few per cent, as demonstrated in Fig. 1. 



C. Setup for the numerical simulations 



Using standard Wilson discretization and denoting lattice quantities with hats, the lattice 
action corresponding to our effective theory reads 



Sa = Sw + Sz + V{^,U), 

1 



x,i<j 



(C.l) 
(C.2) 



Sz = 2 X^Tr [n' - Hx)Ui{x)fi{x + i)ul{x) 

^ x,i 

^ / rr. A 



V 



where (3 is the lattice coupling constant 



/9 



agl 



(C.3) 
(C.4) 

(C.5) 



corresponding to a lattice spacing a. 

The superrenormalizability of the theory allows the matching of the fields and coef- 
ficients of the lattice and continuum theories to be performed exactly up to the desired 
order in a using lattice perturbation theory. To 0{a^), which is needed for performing 
simulations with physical MS parameters, the lattice-continuum relations read 



S = 53^^, n = 5311, a 



(C.6) 



and 



^1/93 - ^(6ci +3C3)^2 



+ 



1 



+ 



167r2 
1 

167r2 



S 

gla 

(48c? + 12ci - I2C3) 
(80c^ + 4ci - 4OC2) 



6 








log — + C 


-12C3 


[¥-] 


} 


«/"3d 









+ 0{a) (C.7) 











log^ 


K 


-40c2 


[¥-] 


«/W3d 









+ -ttT. - 16(5 + p)+ 8ki - 4k4 
2 3 



0{a) 



(C.8) 
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where Eq. (|C.6| ) follows from tree level scaling and Eqs. ( C.7 ) and ( |C.§| ) from two- loop 
calculations, and in the latter two the numerical coefficient S has no relation to the cor- 
responding field in our effective theory. For details of a corresponding calculation with 
the gauge group SU(3) and the definition of the lattice constants, see Refs. [8] and [21]. 
Finally, plugging in the numerical constants and setting /isj = 53, we get 

f , , 4 2.38193365 , ^ . , ^ 
^1 = ^1/53 (2ci + h)P 

+ {(48c? + 12ci - I2C3) [log 1.5/3 + 0.08849] - 6.9537 C3} + 0(/3"^)), (C.9) 
n 7Q3Q77Q 

h2 = b2/gl - "-^'f '^^ 10C2 + C3 + 2)/3 (C.IO) 
47r 

+ Y^|(80c| +4c| -4OC2) [log 1.5^ + 0.08849] - 23.17895 C2 - 8.66687| +0{p-^). 
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